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We present exact and asymptotic results for clusters in the one- dimensional totally asymmetric 
' exclusion process (TASEP) with two different dynamics. The expected length of the largest cluster 

is shown to diverge logarithmically with increasing system size for ordinary TASEP dynamics and as 
. a logarithm divided by a double logarithm for generalized dynamics, where the hopping probability 

— ' of a particle depends on the size of the cluster it belongs to. The connection with the asymptotic 

theory of extreme order statistics is discussed in detail. We also consider a related model of interface 
' growth, where the deposited particles are allowed to relax to the local gravitational minimum. 

(N 

Q . I. INTRODUCTION 

The totally antisymmetric exclusion process (TASEP) is one of the most studied models of nonequilibrium phe- 
nomena and has attracted both physicists and mathematicians. Many of its features, including density profiles and 
currents, have been solved with periodic and open boundary conditions In particular, the so-called matrix 

formulation has provided an elegant way of exploring the properties of the steady state and has been successfully 
applied to get nontrivial results for the so-called shock solutions and for systems with several species of particles . 
Recently, it has been used to solve the exact large deviation functional for the density profile in the case of open 
boundaries [||. 

TASEP has been used to model free traffic flow Q| and traffic jams induced by a crossing j8|. Traffic interpretation 
naturally raises a question of queue lengths in the system. Although this problem is equivalent to one-dimensional 
site percolation in the canonical ensemble, if the steady state with usual exclusion dynamics and periodic boundary 
conditions is considered, being exactly solvable, it may still give some insight to problems, which are too complicated 
to solve in more complex situations. The statistics of the longest queue, or particle cluster as we shall call them, 
I ' is interesting because it measures the size of the large deviations in the system. Largest clusters in percolation in 
general dimension were recently studied by Bazant |J. An exact solution, with limiting behaviors, is presented in this 
• article for the one-dimensional exclusion process. In exclusion processes, one has also the possibility of changing the 
t-H | microscopic dynamics, which leads to percolation problems with different types of cluster statistics. For example, a 
model for 'reckless drivers', who are trying to escape from their chasers is considered in the present article. We show 
that this dynamics results in Poisson statistics for the cluster sizes, which in turn makes the properties of the largest 
cluster very different from those in the case of the usual TASEP dynamics. 

TASEP is also related to two other well-known nonequilibrium systems, namely interface growth and directed 
polymers in random media 0. In particular, the cluster properties are translated to the excursion properties of the 
' growing interface and the largest clusters measure the long excursions before turning back. 

This article is organized as follows. In Sec. ||, we present the basic definitions and the mapping between TASEP 
and the so-called zero-range process. The properties of clusters in stationary TASEP are studied usin g thi s ma pping 
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in Sec. III. In particular, we present analytical results for the largest particle and hole clusters in Sees. 
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In Sec. IV extremal statistics of clusters for TASEP with the modified dynamics is considered. Finally, in Sec. 



discuss a related model of interface growth. 
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II. THE MODEL 



X 

The model system we analyze in Sees. [II and IV is the totally asymmetric exclusion process (TASEP) with n 



particles hopping to the right in a one-dimensional lattice of N sites, where N is an even number, with periodic 
boundary conditions. For convenience (see below), we first restrict 

p := n/N = 1 - 1/k, where k=2,3,... . (2.1) 

The dynamics, however, differs from that of the ordinary TASEP. We introduce another parameter z — 0,1, ... ,n, 
which for z > 1 makes the hopping probability to depend on the length of the queue behind the particle in question 
as follows: After randomly picking an occupied site, check how many occupied sites there are in between the chosen 
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site and the next empty site. If this number is less than or equal to z, the rightmost particle in the cluster makes a 
move to the right - otherwise nothing happens. 

Now the case z — is clearly the usual TASEP. In the present article, we concentrate mainly on the cluster length 
properties of the two cases z = and z = n. We shall also comment briefly on a model of interface growth, where the 
deposited particles always relax to the local minimum, which is equivalent to a certain symmetrization of the z = n 
case. 



A. TASEP and the zero-range process 



Many properties of the exclusion process can be accessed by using the mapping [|10 11 1 to the zero-range process 
introduced in Ref. |T^| . The mass variable M,, i = 1, 2, ...N — n, of the zero-range process is the number of particles 
between i th and (i + l) th hole in TASEP. Since in TASEP there are n particles and N — n holes on a ring of N sites, 

there will be N — n mass variables each attaining integer values from zero to n with the constraint that 53 j=1 n Mi = n. 
A jump of a particle in TASEP corresponds to a change (Mj,M,+i) — > (Mj — 1,-Mj+i + 1) for some i. This way 
TASEP is reduced to the problem with n particles hopping to the right on a lattice of N — n sites, now with multiple 
occupancy allowed. 

This mapping is, however, not one-to-one, because there are N/(N — n) TASEP configuration corresponding to one 
configuration of the mass variables. For example, the number of all mass variable configurations equals the number 
of solutions of the equation Mi + M<x + . . . + Mn- h = n in non-negative integers. This number is (mJZ-i)- But 
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n(jv- -i) = ( )' w hi cn equals the number of all corresponding TASEP configurations, i.e. the number of random 
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walk paths from (0, 0) to (N, N — 2n) . Note that without the restriction of Eq. ( |2.l| ) this relation between the numbers 
of configurations would have applied on average only. 



As discussed in Ref. |10|, the joint stationary distribution of the zero-range process corresponding to TASEP with 
general z is the product of the marginal measures of the individual sites. For the cases z = and z = n, however, the 
stationary state is particularly simple and can be constructed in a straightforward way. 



III. SIMPLE EXCLUSION PROCESS, z = 



A. Stationary cluster size 

It is general knowledge that the stationary distribution of TASEP on a ring of N sites is the uniform distribution [Q , 

the probability of each distinct configuration being ( ) . This follows immediately from the fact that the transition 
probability to a given state equals the transition probability from the state in question to other states, which makes 
the transition matrix doubly stochastic |l3| . The mapping from TASEP to the zero-range process being one-to-one 
up to a multiplicative constant, the stationary distribution of the zero-range process corresponding to z — is also 
uniform. 

Due to the translational invariance of our model system the mass variables Mi are identically distributed and the 
general behavior is described by local expectations. The stationary distribution for size Xi of the i th cluster in TASEP 
can then be expressed in terms of the stationary mass variable distribution by 

P(X t =k) = P(Mi - k\Mt > 1). (3.1) 



Next we use translational invariance and the uniformity of measures to conclude that Eq. (3.1) equals the number of 
solutions to 

Mi + Af 2 + . . . + Mtf-n-i =n-k, (3.2) 
where Mj's are non-negative integers and k > 1, divided by the difference of f^ -1 ) and the number of solutions to 



Eq. (3.2) with k = 0. Therefore 
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n 



Note that one could have derived this result also by simply counting the random walk paths starting with exactly k 
steps up and ending with a down step and normalizing this number by the number of random walk paths with the 
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first step up and the last step down. From Eq. (B.3) we immediately obtain the limiting distribution for the cluster 
lengths, 
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P oc {X l = k) = 

in agreement with the fact jjj that in an infinite system size the measure for TASEP is the product measure. 
The expected cluster size is 



(3.4) 
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and its variance is 



VarX = 



(p-l/N)(l-p-l/N)(l-l/N) 



(l-p) 2 (l 
P 



(\- P r 



-p + l/N) 
as N 



oc. 



(3.6) 



The remarkably simple result of Eq. (3__5) shows that the expected cluster size converges to a finite value and it 
is the number of clusters that has to diverge in the N — > oo limit. This can be also verified by the following direct 
evaluation. Let us first calculate the probability that there are k clusters in the system, i.e. that there are k occupied 
sites in the mass picture. If the first site is to be occupied, place k sites on a ring, label one of them 1 and fill them 
with n identical objects, such that every site gets at least one object, in ways. After that, place N — n — k 

vacant sites in between the k occupied sites in (^jT"^ 1 ) ways. On the other hand, if the site 1 is to be vacant, place 
k + 1 sites on a ring, label one those 1 and fill the sites 2, 3, . . . , k, such that every site gets at least one object, in 



(fell) WS- Then add N 
distribution, 



k — 1 vacant sites in between in ( 



ways. Thus we get the hypergeometric 



P(#C = k) 



which immediately yields 
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From Eqs. (3.5) and (3.8) we find the intuitively expected relation (#C) ■ (JQ 



n. 



B. Length of the longest cluster 

The knowledge of expectations of maximal objects is of great importance since it measures the size of large deviations 
in the system. In this section we study the statistics of the longest cluster, max^ JQ, in TASEP. In particular, we 
discuss the finite size and limiting distribution functions, related expectations and the effects due to discreteness of 
the sample space of the mass variables. 

The cumulative distribution function (c.d.f.) 



F PtN (k) := ?(maxl, < k) = P{ max Mi < k) 



Ki<N-r 



(3.9) 



can be calculated as the ratio of the number of solutions to 



Mi + M 2 



M 
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(3.10) 
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in non-negative integers with Mj < k and the number of its solutions with Mj < n, Vi, respectively. The first of these 
two numbers equals the coefficient of the n th order term in 



(l + a + ...+s*)"-"=(±Z£_) (3.11) 



1 - x k+1 



N-n 



and therefore 



where 

l*nM k ) := min {^ ~ n, [n/(k + 1)J} (3.13) 
and |_-J stands for the integer part. Note that for k > n/2 the distribution function takes the simple form 

\N-n-l ) 

(V) 

The expectation calculated in terms of the c.d.f. is 



F p!N (k) = l-(N-n)^^. (3.14) 



fc=i 

so that, after simplifications, it can be written as 



(maxl,) = V k ■ P({Mi <k\fi}n {3j : Mj = k}) 
k=l 

n 

= ^2 [l - P(Mi < k - 1 Vi) , (3.15) 



<m ^ i>= cvjS 5 (_ ) ( ^ )(iv-n-i)' (3 - i6) 



where /i n ,jv(fc) was defined in Eq. (|3T3|). 



C. The limiting extremal distribution 

It seems difficult to extract the limiting distribution and the type of the divergence of the mean from the exact 



results of Eqs. (3.12) and (3.16) directly. For that purpose, we next consider an alternative approach based on 
independence. 

Intuitively one would expect that as the system gets larger the weaker become the correlations between the clusters. 
In fact, it is easy to see that the mass variables are asymptotically independent, in that, for k fixed, 

k 

P(M n =oi,M ia =a 2 ,...,M lk =a k ) — ►JJP^Mj, = o,), (3.17) 

z=i 

where i; ^ i m , V7,m, and N — > oo. Therefore, we can approximate the c.d.f. F Pi N(k) by using the product measure, 
i.e. 

P(maxXi < A;) ~ P(M 1 < k) N - n . (3.18) 

i 

By writing 
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P(Mi <k) = 1 - 

— ► 1 - p k+1 as N -> oo, (3.19) 
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we get in the large N limit with p fixed 
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-(fc+l)logl/p 



1 - e 

fclogl/p + logNp(l-p) 
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(3.20) 



This result is, for fc sufficiently large, in agreement with the earlier result for the tail because Eq. (3.14) yields in the 
large system size limit 



F p , N {k)~l-N{l- p)e-^ io ^, 



(3.21) 



which is simply the first order approximation for Eq. (3.20) 



The c.d.f. in Eq. (3.20) can now be expressed in terms of a scaled variable 



logjVp(l-p) 
log 1/p 



(logl/p)" 



(3.22) 



From the asymptotic theory of the extremes for independent and identically distributed (i.i.d.) variables [[14 15 we 
then obtain the Gumbel, or Fisher-Tippett, distribution 



P (Z p ,n < y) ~ e 



(3.23) 



Strictly speaking, however, the function on the right hand side is continuous, whereas the c.d.f. on the 
left must be piecewise constant with an intrinsic thi cknes s log 1/p, so that the relation holds only when 
(y + log Np(l — p)) /log (1/p) is an integer. In fact, Eq. ( 3.23 ) would be the correct result for the continuous ex- 
ponential parent distribution, whereas it has been shown by Anderson ]Tr| that the discrete geom etric distribution 
does not belong to the domain of attraction of the Gumbel distribution at all, and that Eq. ( 3.23 ) should indeed be 
replaced by 



limsupP(Z p jv < y) < e 

N^oo 

liminf P (Z Pt N < y) > e~ 



-y+log i/p 



(3.24) 
(3.25) 



i.e. the limiting distribution function has Gumbel envelopes. 

The expectation of a random variable with the Gumbel distribution is known to equal Euler's constant 7 rs 0.5772 
fill , but the evaluation of the expectation of a discrete variable with Gumbel envelopes is not easily found in the 
literature. A direct calculation using the proper point probability function seems tedious and, therefore, we argue on 
heuristic grounds that since the piecewise constant c.d.f. of the scaled variable Z p ^ is confined between two envelope 
functions with expectations 7 and 7 + log 1/p, the final result should be written as 



(max X;} 



7 + logXp(l - p) 1 



log 1/p 



2' 



(3.26) 



where N is large. Fig. |l|(a) shows that the exact result Eq. ( [3.16 ) and the asymptotic formula Eq. (3.26) are in 
excellent agreement with the simulation data |p7| . 

The effect of discreteness of the mass variables can be seen even more clearly in the behavior of the variance 
than of the expectation. For continuous exponential variables, the variance of the maximum converges to a constant 
<7 2 7r 2 /6 [ pi , where the scaling parameter a equals (logl/p) -1 in our case. The simulation results of Fig. |l|(b) display 
persistent fluctuations around this value, with the distance between the maxima diverging logarithmically in N. 
This kind of behavior has recently been observed also by B azant || in connection with percolation theory and the 
limiting distribution functions described by Eqs. ( [3.24 ) and (3.25) were called Fisher-Tippett limit cycles because of 
the quasiperiodic fluctuations of the variance. 



IV. TASEP WITH ASYMMETRIC DYNAMICS, 



We now turn to the case z = n. This means that every time we pick a particle at random, the rightmost particle in 
the cluster, where the chosen particle belongs to, makes a move. In the mass picture this translates to n independent 
random walkers on a ring of N — n sites pTo|] . 
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A. Stationary cluster size 



The joint stationary distribution for n independent random walkers on a ring of N — n sites is 

P(Mi = mi, . . . , Afjv-„ = mw-n) = 7T7 1 ( " ) , (4.1) 



(iV — n) rl \mi . . . rriM-i 

where Mi is the number of walkers on the site i. By Eq. ( |3.l| ), the cluster length distribution for the TASEP variables 
Xi is now 

pN-k 



P(X l = k) = ± ^ , (4.2) 



1 



(* iV(l-p)) 



^ v fc / ' 



the limit of which is Poisson with the proper normalization, 



\ — k 

P QO {X l = k)= where A = (4.3) 

re! (e A — 1) 1 — p 

The parameter A introduced above describes the length scale of correlations in the system. Remember that in the 
case z = the corresponding parameter was (logl/p) -1 > A, so that the correlations are expected to die out faster 
in the present case. Now we have the mean and the variance of the length of particle clusters 

(Xi) = - 



pN 

1-(1- 



N(l-p) 

A 



1 - e 

VarX, = (Xi) (1 + A - (X t )) 
A 



—A 



as N -> oo, (4.4) 



1 + A r as N -> oo. (4.5) 

1 — e A / 

Since there are ( N ^ n ) ways to choose k sites from a ring of N — n sites and k\S(n, k), where S(n, k) is the Stirling 
number of the second kind, ways to put n particles into k sites such that each site gets at least one particle, the 
distribution for the number of clusters now reads 

which again gives (#C) • (Xi) = n. 

Due to the lack of particle-hole symmetry, the properties of the hole clusters differ from those presented above. 
However, the distribution for the length Hi of the i th hole cluster can be easily constructed in terms of the distribution 
for the mass variables, 

, , f P(M i+l > 0) iffc = l 

^ [Tlt \ P(M i+ i = 0, M l+2 = 0, . . . , Mi+u^x = 0, M i+k > 0) if k > 2 

Jfe-1 \ pN / fc ^ w 



The corresponding c.d.f is 



-> (e A - 1) e~ xk as iV oo. 



P(Hi < k) = 1 - I 1 - — ) — > 1 - e~ Afe as iV -> oo, (4.8) 



A(l - p) 



and the expected length of hole clusters 



N—n , h —\ \ pN 1 



fe=i 



G 



B. Length of the longest cluster 



The exact longest particle cluster c.d.f. for the case z — n reads 



Pfmaxl, < k) = , Ar 1 . V ( " ), (4.10) 

s>«=- 



which seems to be, from the viewpoint of applications, even less instructive than Eq. ( 3.12 ) was f or the ordinary 
TASEP. However, an approximate form for the tail of the distribution can be obtained from Eq. (4.10) assuming that, 
for k sufficiently large, at most one of the masses exceeds k, 



3j: 

(iV-n) V ( " Y N ~"r \n (4.11) 

1 ' ^ \m N -J (N-n) n 1 ' 



x £ 

mi <fc. 



n - mAr-7 



mi ... mjv-n-i/ (-/V - - l) n_rriw - 



Here the sum in the last term is the probability that on a ring of N — n — 1 sites with a total of n — m^-n independent 
random walkers there are at most k walkers on a single site. But for k > (n — l)/2, and approximately perhaps even 
for smaller k values, this probability equals unity and therefore we have 



PfmaxXj < k) » 1 - V 



n \(iV-n-l) 7 



TOjv-«y (N-n) n 1 



1 - iV(l - p)A fc e- A V — ^ - 

\ k er x 

^"wnrv (4 - 12) 



when k is sufficiently large and also N large at the last two stages. 



Again, one can show that the cluster sizes are asymptotically independent, in that Eq. ( 3.17 ) holds with 

W = k) = (N (N-nr~ k (fc) -> PM = k) = T aS N ^ °°' (4 ' 13) 



Therefore, Eq. (4.12) suggests that the limiting distribution function is, as in the z = case, approximately given by 
the product 




(4.14) 
(4.15) 

where k > A and TV is large. From the mathematical point of view, however, one should be careful with this 
approximation, because the distribution function for the maximum of i.i.d. Poisson variables does not actually converge 
to the Gumbel distribution - not even in sense of continuous Gumbel envelopes as it was in the z — case p6| , p^| . 
Namely, Anderson has shown Jl6[ that the probability of the maximum concentrates on two consecutive integers, i.e. 
there exists a sequence of integers /„ (A) such that 
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lim P( max M l = I n (X) or I n (X) + 1) = 1. 

n/X — >oo l<?<n/A 

The sequence (/„) then obviously determines the type of divergence of the mean. It was shown by Kimber 
in the leading order this sequence diverges as 



/n(A) 



logn 
log log n ' 



(4.16) 
| that 

(4.17) 



However, we would like to remark that, despite of its shortcomings, the approximation of Eq. (4.15) yields the same 
functional form for the expected length of the longest cluster: In the leading order 



(max Xi 



n-l 

^ t X{Km x } i 
k=0 



(4.18) 



where X{ } is the indicator function and m x is some crossover mass, where the distribution function increases most 
rapidly. Eq. (4.15) then gives the very interesting result 



(m x + 1)! 

\m x 



e~ x n, 



(4.19) 



i.e. the expected length of the longest particle cluster diverges as the inverse relation of Eq. (4.19) in n. In the special 
case A = 1, i.e. p = 1/2, the type of divergence reduces to inverse factorial. Furthermore, by Stirling's formula one 
obtains the same functional form as in Kimber's result in Eq. (4.17). Our analytical result is again supported by 
simulations as seen in Fig. ||(a). 

The plot of the variance in Fig. ^(b) sho ws aga in fluctuations, and the divergence of the distance between maxima is 
now dictated by the inverse relation of Eq. ( 4.19 ). In general, it can be seen that the fluctuations are more pronounced 
than those in the ordinary TASEP. We also observe that the variance has a global maximum at a finite system size. 
Below this maximum, the correlation length A is comparable to the system size, which results in large fluctuations of 
the maximum. For N 3> A the system can be considered to consist of a large number of independent copies and the 
i.i.d. beha vior is recovered. In the end, one should notice that, unfortunately, the continuum approximation of the 
form Eq. ( 4.1 5| ) with the factorial replaced by the gamma function cannot be used to evaluate the variance: In the 
continuous case, the maximum concentrates on one real number, in that the properly scaled variable converges to a 
distribution degenerate at zero |T^ , |l6|| . 

In Fig. |^ we show the expectation and the variance also for the hole clusters. According to Eq. d4.7| ), the limiting 
distribution for the lengths of the hole clusters is again geometric, as in the case of z — dynamics, but the number 
of clusters is random and the length of the longest cluster depends on that number. The solution to this pr oblem can 
be found by defining new mass variables M* as the number of holes between i th and (i + l) th particle. Eq. (4.8) then 
tells that 



P(M* <k) = l-C, 



p,N 



1 - 



N(l - p) 



pN 



C p e 



-Xk 



as N — > oo, 



(4.20) 



where (C P] at) is a sequence of constants with limC Pj Ar = C p > 0. Since there are now n new mass variables and the 
parent distribution is geometric-like, the limiting extremal distribution has again Gumbel envelopes and the divergence 
of the mean is logarithmic, as seen in Fig. ||(a). 



V. DISCUSSION 



We shall now discuss our results from a different point of view. There is a well-known mapping from TASEP to 
a certain lattice model, namely the one-dimensional single-step model, of interface growth In this mapping each 
particle is considered a unit step down and each vacancy a unit step up. The function h(x,t), where x is the spatial 
and t the time coordinate, obtained this way is defined to be constant between the integers and, say, right-continuous. 
Clearly, the condition p ^ 1/2 implies existence of a global tilt. The lengths of the clusters in the TASEP picture 
translate to the lengths of the decreasing parts of the interface. Note that the z = case has particle-hole symmetry 
and the results obtained for the decreasing parts (corresponding to particle clusters) are valid for the increasing parts 
(hole clusters) as well. For z > 1 the properties of the increasing parts must be calculated from the hole distribution. 

The z > 1 dynamics can be seen to be equivalent to deposition rule such that the deposited particles relax towards 
the local minima from the decreasing parts of the height function h(x, t). If the particle lands on an increasing part, 
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nothing happens. Since this kind of asymmetry is not very common in real applications of the interface interpretation, 
we also studied numerically a spatial symmetrization of the z — n case with p — 1/2. This symmetrization belongs to 
the Edwards- Wilkinson universality class |Q. In it, deposited particles always flow downhill to a local minimum. In 
case the particle lands on a hill top, the minimum is determined using a fair coin. The results are plotted in Fig, |g . 
The mean length of the longest cluster diverges now faster than in the asymmetric z = n case studied in Sec. IV, 
the earlier analytical result for which is shown for comparison in Fig. ^(a). This is consistent with the fact that on 
average the hole clusters are considerably shorter than in the case of asymmetric dynamics, which is compensated 
by longer particle clusters. The other features can be seen to be quite similar to those of the asymmetric case. For 
example, the variance of the longest cluster has a maximum at finite N and quasiperiodic oscillations as before. 

In conclusion, we presented analytical and numerical results for the properties of the clusters in asymmetric exclusion 
processes with two different dynamics. In particular, we studied the properties of the longest cluster in the system 
and showed that the large system behavior agrees with the asymptotic theory of the extremes for independent and 
identically distributed variables. The expected length of the longest cluster was found to diverge logarithmically with 
increasing system size for the ordinary TASEP and as logarithm divided by double logarithm in the case of modified 
dynamics, which corresponds to independent random walkers in the mass picture. In the latter case, the length of 
the longest cluster is distributed among two consecutive integers in the large system limit. 
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FIG. 1. (a) Expectation value of the size of the longest cluster and (b) its variance as a function of the system size N fo r 
z = with p = 1/2. In (a), the dotted curve is the exact result Eq. (3.16) and the full curve the asymptotic formula Eq. (3.26). 
Crosses denote the simulation data. 



10 




FIG. 2. (a) Expectation values of the size of the longest particle cluster and the size of the longest hole cluster and (b) 
their variances as a function of the system size N for z = n with p = 1/2. Crosses and squares denote the simul ation da ta fo r 
particle clusters and hole clusters, respectively. In (a), the dotted curve was obtained by substituting Eq. ( 1.14) in Eq. (3.15) 
and the full curve by taking the crossover mass to be the point where the approximate c.d.f. of Eq. (4.15) obtains the value 
1/2. The dash-dotted line is the function 0.987 log iV, where the prefactor was chosen to produce the 'best fit'. 
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FIG. 3. (a) Expectation values of the size of the longest cluster and (b) its variance as a function of the system size for 
the symmetrized model with z — n and p — 1/2. Crosses denote the simulation data for clusters in the symmetrized model. In 
(a), the asymptotic curves for the corresponding asymmetric case from Fig. ^ are shown for comparison by the dashed lines. 
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